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Abstract 

Axion production due to photon-axion mixing in tangled magnetic field(s) prior to re¬ 
combination epoch and magnetic field damping can generate cosmic microwave background 
(CMB) spectral distortions. In particular, contribution of both processes to CMB p distor¬ 
tion in the case of resonant photon-axion mixing is studied. Assuming that magnetic field 
power spectrum is approximated by a power law Ps(fc) oc k n with spectral index n, it is 
shown that for magnetic field cut-off scales 172.5 pc < As < 4 x 10 3 pc, axion contribution 
to CMB p distortion is subdominant in comparison with magnetic field damping in the cos¬ 
mological plasma. Using COBE upper limit on p and for magnetic field scale A b — 415 pc, 
weaker limit in comparison with other studies on the magnetic field strength (Bo < 8.5 x 10 -8 
G) up to a factor 10 for the DFSZ axion model and axion mass m a > 2.6 x 10 -6 eV is found. 
A forecast for the expected sensitivity of PIXIE/PRISM on p is also presented. 
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Introduction. During the last decades have been done intensive studies regarding the exis¬ 
tence and nature of primordial magnetic field(s) at both small and large scales. Its existence 
could have strong impact in different cosmological scenarios such as bing bang nucleosynthesis 
(BBN), structure formation, CMB temperature anisotropy etc. In general, in all those scenar¬ 
ios, it is possible to probe its existence only indirectly, namely through the coupling of magnetic 
field with the cosmological plasma. Consequently, based on information that we have on BBN, 
CMB temperature anisotropy etc., it is possible to speculate about the magnetic field structure 
and estimate its strength at a given scale. In particular, CMB temperature anisotropy has been 
one of the most important benchmark to test the existence of primordial magnetic field(s). In¬ 
deed, an ubiquitous, anisotropic and homogeneous magnetic field with strength at present time 
B 0 < 3 X 10 9 G would create the observed CMB temperature anisotropy due to anisotropic 
expansion of the Universe [1]. For a general review on cosmological magnetic field see Refs. [2j. 

In the presence of a large scale magnetic field, CMB photons can in principle convert into 
axions or other similar particles due to their coupling with the magnetic field. In Refs. |3j we 
have studied such a mechanism in the presence of large scale uniform (spatially homogeneous) 
magnetic field and applied it to CMB spectral distortions and temperature anisotropy. However, 
several interesting questions arises as to what happens in the case when the background magnetic 
field is not homogeneous (tangled magnetic field). Does the magnetic field has an impact on the 
spectral distortions? Is the impact of the magnetic field on the spectral distortions dominant 
or subdominant with respect to photon-axion oscillations? 
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As in the case of density perturbations in the primordial baryonic plasma suffering from Silk 
damping, we can expect that a spatially varying magnetic field can couple to the baryon plasma 
and dissipate energy. This would eventually lead to the damping of the primordial magnetic 
field spectrum in different scales [3|. In general, a distinguishing feature of non homogeneous 
magnetic field in comparison with an uniform field is that the former can have an impact on the 
CMB by distorting its spectrum. Indeed, it has been shown in Ref. [5] that a spatially varying 
stochastic magnetic field may significantly dissipate in the cosmological plasma prior to recom¬ 
bination epoch. By dissipation, the magnetic held energy would transform into kinetic energy 
of cosmological plasma and in turn plasma’s kinetic energy would be efficiently transformed 
into heat due to high shear viscosity of plasma. In the limit when photon mean free path Z 7 
is smaller than magnetic held mode A, -C A, Alfven, slow and fast magnetosonic waves with 
A < d 7 are effectively dissipated where d 7 = (fit) 1 / 2 is the photon diffusion length and t is the 
cosmological time. 

In the case when there is an energy injection into the cosmological plasma such as con¬ 
version of magnetic held energy into heat, electrons gain energy and the electron temperature 
T e becomes higher than the photon temperature T, T e > T. Depending at which redshift the 
magnetic held energy is converted into heat, this effect would eventually lead to CMB spectral 
distortions if energy is injected for redshift z<2x 10 6 . For an early treatment of CMB spectral 
distortion see Ref. [6j, for further developments see Ref. [3 and for production mechanisms of 
spectral distortions see Ref. [8]. 

Dissipation of tangled magnetic field. During the evolution of the universe, it is usually 
assumed that the conductivity of the cosmological plasma is infinite. In this case the held 
amplitude scales as B ~ Boa~ 2 (t) where a is the cosmological scale factor. Even though 
this is a good approximation, it does not reflect the more general case, namely in the case of 
tangled magnetic helds when the magnetic held dissipate energy. In order to make contact with 
our results that follow, we assume that the magnetic held is generated by random processes 
(stochastic) in the early universe during inflation or radiation epoch and it evolves according to 
the following law 

B(x,i) = b 0 (x,i) , (1) 

where ao is the cosmological scale factor at present epoch and bo(x, t) is the tangled mag¬ 
netic held of which amplitude evolves as b = bo exp(i dt'u) where ui is the magnetic held 
mode frequency, see Refs. [3j and [5|. Moreover we assume that magnetic held is statistically 
homogeneous and isotropic with ensemble average 

(bfi k)6*(q)) = J 3 (k - q)i^(k)P B (A;), (2) 

where P,j is a projection tensor and Pb is the power spectrum of primordial magnetic held 
that in general is assumed to be a power law, Pb = Ck n with C a constant and n the spectral 
index of the magnetic held. The constant C is hxed by taking the spatial average (or ensemble 
average) of the energy density of the magnetic held over a volume V 

MM = (B ° (X) 2 B0(X)> = f, (3) 

with comoving cut-off wavelength A b = 27T /k\ where k\ is a cut-off comoving wave-vector. The 
cut-off wavelength (or wave-vector) is in general a free parameter that is connected with Fourier 
decomposition of B(x) and in principle can assume values from zero to infinity. However, for 
physical reasons it depends essentially on the generation mechanism of the primordial magnetic 
held. In the case of magnetic held generated by casual mechanism(s), the value of A# should 
be smaller or equal to the Hubble horizon while in the case of magnetic held generated by non 
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casual mechanism(s) (in general negative spectral indexes) the value of A# can be greater than 
Hubble horizon. As we see below, in the case of casual mechanisms, we consider that As to be 
smaller or equal to Hubble horizon during the p epoch, see Figs. 01 and [3j In the case of 
magnetic fields generated by non casual mechanisms we set for simplicity Ag to be of the order 
of Mpc even though it can be larger than this value, see Figs. |4b| and [5j 

With this kind of normalisation the magnetic field power spectrum is given by 

Another possibility on fixing the constant C, that is also used in the literature, is to use a 
Gaussian filter e~^ k ^ kx ^~ in the definition of pg(io)> see Refs. [9]. In this case the form of Pg is 
different from Eq. Q but the spatial average of pg(to) remains invariant as it should be. In 
this work we shall not adopt this definition. 

Let Qb indicate the energy loss per unit time of magnetic field that would convert into heat in 
the plasma. If energy injection (or heat) occurs in the redshift interval 2.88 x 10 5 < z <2x 10 6 , 
the Compton scattering would eventually create a Bose-Einstein distribution for the photon 
spectrum with chemical potential /Q The evolution of chemical potential with respect to time 
is governed by the Sunyaev-Zel’dovich equation | 6 ] 


dp = l±_ 1 1 Qb 

dt tdC PR ’ 


(5) 


where tdc = 2.09 x 10 33 (1 — Y p /2) 1 (/; 2 Hg) 1 (l + z) 9 / 2 s is the characteristic time for double 
Compton scattering. Here Y p ~ 0.24 is the helium primordial weight by mass and /i 2 Hg ~ 0.022 
is the density parameter of baryons {TO]. 

To solve Eq. ([5]) we need to know the rate of heat flow into the plasma due to magnetic 
field dissipation. It can be shown that in the photon diffusion limit i.e. T ^ 1 <C Ag [3| (where 
r 7 = a T n e with eg being the Thomson cross section and n e the number density of free electrons) 


Qb = Bl (n + 3) f kx k n+A f 2k 2 1 A 

PR 2 p R {to) k™ +3 Jo 5(1 + z)r 7 (f 0 ) 6XP \ k 2 D (t 0 ) {l + z) 3 )' 


where pr is the energy density of relativistic particles, to denotes the present time and k 2 D (to) = 
15r 7 (t 0 )/G — 6.27 x 10~ 19 r 7 (t 0 ) s^ 1 . The term f* is connected to the thermalization redshift^J 
z fl , through the relation t* = htdc/QQ + 2 ^) 5//2 ). Substituting Eq. ([ 6 ]) into Eq. ([5]), the general 
solution of Eq. ([5]) in terms of the redshift is given by 


p(z) = 


1.4 P 2 


(n + 3) 


10 / 9g(to)r 7 (t 0 )^ +3 io 


rk a rz 

Jo Jz 


dz'dk 


k n+4 

(1+P) 4 


exp 


1 + z' 
1 + z u 


5/2 


exp 


2 k 2 


k 2 D {t 0 ){l + z’) 3 


( 7 ) 

where Z{ is an initial redshift, Zi S> z p and the term proportional to p(zi) is absent since for Zi 
we have p(zi) = 0. In obtaining Eq. ([7]) we have used the fact that in the radiation dominated 
universe dt = — dz 2t* / (1 + z) 3 and z is the redshift in the radiation dominated universe z <C z p . 

In general is not possible to find analytic solution of Eq. ([7]) due to the non trivial form of 
the integrands. Indeed, one can recognise that the double integral of Eq. 0 can be expressed 
in terms of the incomplete gamma functions, 7 (s, 2 A: 2 /A’| ) ) where s is an integer that in our case 


1 The chemical potential introduced here is an dimensionless quantity and is related to the thermodynamical 

chemical potential pther by p = —pther/T. 

“The thermalization redshift is the redshift that for z > 27 the CMB spectrum is in thermal equilibrium and 

for z < the spectrum is a Bose-Einstein distribution, see Refs. J 2 - 
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is either s = (n + 5)/2 or s = (3n + 9)/5. However, it is possible to consider some limiting cases 
that allows to find analytic expressions in terms of Euler gamma function E. Let us consider 
the limit k\ k 2 D (to)(\ + z M ) 3 and then evaluate the residual chemical potential at redshift 
z = 0 (today). In this limit we get the following relation between magnetic field strength Bq 
and p 

nr / kn \ ~( n + 3 )/ 2 

B °= 3 - 19 x i° v^(e) g - (8) 

where C n is a constant 

C n = 1.4 r(n/2 + 5/2) T(3n/5 + 9/5) 2" (n+5)/2 (6/5) (n + 3), (9) 

Here the term &D = kr)(to)zf/ 2 in Eq. ([8]) is the scale damped by one e-fold at redshift z /t . Its 
corresponding co-moving wavelength is \d = 2ir/kD = 415.5 pc. On the opposite, in the limit 
k\ < k 2 D (t 0 )( 1 + z^) 3 , we get 

B„ = 3.19 X (^) HO) 

where D n is a numerical constant that is given by 

^ = 1.4T(-6/5) (^|) (6/5). (11) 

Axion contribution to spectral distortion. We have seen that tangled magnetic fields can 
dissipate energy and create p distortion in the early universe. However, their presence would 
make possible the transition of CMB photons into axion^J In Refs. |3] we have derived the 
equations of motions for the photon-axion system in the steady state approximation in the case 
of uniform magnetic field. Here we calculate the transition probability in the resonance case in 
presence of tangled magnetic field. The resonant regime is the axion mass range that makes 
resonant transition in the redshift interval 2.88 x 10° <l + z<2x 10 6 , see Ref. [3] for details. 
The difference in this case is that the magnetic field depends on the position, B(x, t). However, 
one does not need to calculate the equation of motion for the density operator p again. It is 
only sufficient to replace in the equations of motions for p , B(i) —> B(x, t) and take the spatial 
average of the transition probability, P a —> (P a ). Since the transition probability depends on 
Bq [3] and using the fact that (Bo(x)Bo(x)) = Bg, in the resonant case we get 

(P a (T)) = 5.75 xW~ 27 xC 2 ai B 2 nG ^ , (12) 

where B n c = (Bq/hG) and C a7 is defined as 




<17 


/E 24 + ro\ l + io 
Viv “ 3 1 + w) w 1 / 2 ’ 


(13) 


where for w = 0.56, \C a ^\ — 4 for E/N = 0 (KSVZ axion model) and |C a7 | 1.49 for 

E/N = 8/3 (DFSZ axion model). Here w is defined in terms of the mass ratio of up and down 
quarks, w = m u /md. For small chemical potential p we can write ( P a ) = p,e x /(e x — 1) where 
x = uj/T with uj being the photon energy and T the CMB temperature. In this case we can 
easily find 

B 0 = 6.76 x 10~ n _ ^ G, (14) 

fTla Cq 7 

3 In this paper we focus only on the QCD axion (hadronic axions) 
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where T and fh a = m a /eV are respectively the resonance temperature and axion mass. Since 
we are looking for spectral distortion in the redshift interval 2.88 x 10 5 < 1 + 2 < 2 x 10 6 , the 
corresponding resonant axion mass is within the interval 2.66 x 10 ~ 6 eV < m a < 4.88 x 10 ~ 5 eV 


m- 

Equation ([ 8 ]) gives only the contribution to /U distortion from magnetic field itself. Now we 
must add to it also the contribution from axion creation from the CMB. Indeed, adding to Eq. 
Q, Eq. (14) we get the following relation between magnetic field strength, ^-parameter and 
Ab 


5 0 = v^f 1.6 x 10- 6 C'- 1 / 2 (Ab/A z? )- ( ^ ) +3.38 x lO- 11 ^^) G > (Ab < A D ) 

\ TUa ^07 / 


(15) 


On the other hand, in the limiting case Ad<Ab and adding to Eq. (10), Eq. (14) we get 


Bo = \fjl ( 1.6 x 10 6 D n , G 2 (Ab/Ad) + 3.38 x 10 11 ————) G, (Ab <C As). 
\ Wla ^<17 J 


(16) 


We notice from Eq. (16) that the magnetic field strength depend on the spectral index n 


only through D n . It is interesting to know at what scales the axion contribution to /x distortion 
is smaller than magnetic field contribution. In the case A b Xd we get 


A B > 4.73 x 10 


while in the case A^ <C Ab we get 


4 1Tl a C a ry 
r»i /2 


2/(n+3) 


A Di (A b Xd) 


(17) 


D 1 / 2 

A b < 2.11 x 10~ J _ n — A d, {Xd Ab). 

Wla'-'a 7 


(18) 


We can see from Eq. and Eq. ( [18] ) that A b does not depend on the average strength of the 
magnetic field Bq but only on n, C a7 and m a . For example for n = —2.9, —2, —1,0,1, 2, 3 we have 
respectively C n = 1.27, 0.78, 0.77,1.1, 2.08,4.93,14.05 and D n = 0.38, 2.71,4.07,4.88, 5.43,5.82,6.11. 
If we consider for example the DFSZ axion model, n = 2 and axions with mass rh a = 3.5 x 10~ 6 
eV we have that for A b ^ 415.5 pc, the axion contribution to /x distortion is subdominant 
to magnetic field damping for Ab > 172.5 pc. In the opposite limit, Ae > 415.5 pc, we get 
Ab < 4.03 x 10 3 pc. On the other hand, if we have /x given by the experiment and ask at what 
scales axion contribution in the Bo — Ab plane is subdominant to magnetic field damping, we 
must simply reverse the inequality signs in both Eq. © and Eq. ( fl8] ) , see Fig. [4aJ 

In Fig. [l]the exclusion plot for the scale averaged magnetic field, Bo vs. Ab is shown. In 
both (a) and (b) the plots for the upper limit on /x found by COBE |11| are shown. Here we 
have chosen magnetic fields with n > 2 which are generated in the early universe by causal 
mechanisms m- For such magnetic fields, the field wavelength A or Ab must be smaller 
than Hubble distance at redshift z, namely Ab < H~ 1 {z). Indeed, in Fig. [I] we have chosen 
H~ 1 (zqcd) < Ab < H^ 1 {z fl ) where H~ 1 {zqcd) ~ 1 pc is the QCD comoving horizon and 
z^ = 2.88 x 10' 5 (lower redshift of /x epoch). Our exclusion and sensitivity plots in Fig. 0 Fig. 

[2] and Fig. [3] have been obtained for /x distortion and n = 2, 3. The region above the solid line 
is excluded with no photon-axion mixing taken into account while regions above the dashed, 
dot dashed and dotted lines are excluded by taking into account it. The exclusion plot with no 
photon-axion mixing has been obtained by using Eq. 
until Ab —> Xd- 


and Eq. (10) and extrapolating them 
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We can see from Fig. 0 Fig. 0 and Fig. [3] that when we take into account photon-axion 
mixing, there are significant deviations for 200 pc < A# < 10 3 pc, in comparison with no photon- 
axion mixing. Depending on the axion mass and axion model, deviations range from a factor 2 
until a factor 11. In our plots we have chosen three representative axion masses, m a = 4.88x ICC 5 
eV, rh a = 1 x ICC 5 eV and fh a = 3.5 x ICC 6 eV. Axions with mass rh a = 4.88 x ICC 5 eV are 
resonantly produced at the beginning of /x-epoch while axions with mass fh a > 3.5 x ICC 6 eV 
are experimentally allowed by ADMX collaboratiorj^] [13]. Axions with mass rh a — 3.5 x ICC 6 
eV give weaker limits on Bq in comparison with axions with bigger masses. This can be seen 
from Eq. (14) where the axion mass is in the denominator. 

In the case of expected limits on /x by future missions such as PIXIE/PRISM [T3], our plots 
are shown as sensitivity plots, see Fig. U Fig. [3] and Fig. [5] For example we can see in Fig. |] 
that PIXIE/PRISM have a much better sensitivity with respect to COBE in the Bq — A b plane. 
Depending on the axion mass and axion model the improvement is in general one or two orders 
of magnitude. In Fig. 4a plots of Eq. (14), Eq. (15) and (16) are shown. In this figure we can 
see the regions where the axion contribution in the Bq — \b plane is dominant or subdominant. 
In Fig. 4b and Fig. [5] plots for negative spectral indexes, for example n = —2.9 and n = —2 are 
shown. In general magnetic fields with negative spectral indexes are generated by non-casual 
processes, for example during inflationary epoch [2J. 



10 100 1000 10 4 


A s (pc) 

(a) 



10 100 1000 10 4 


A s (pc) 

(b) 


Figure 1: Exclusion plot in the parameter space Bq — Xb in the resonant case due to /x-distortion. 
In (a) the exclusion plot for COBE m upper limit on /x and DFSZ axion model for n = 2 is 
shown and in (b) the exclusion plot for COBE upper limit on /x and DFSZ axion model for 
n = 3 is shown. In both (a) and (b) the region above the solid line represents the excluded 
region without photon-axion mixing while the region above the dashed, dot dashed and dotted 
lines represent the exclude region including photon-axion mixing for fh a = 4.88 x ICC 5 eV, 
fh a = 1 x 10 5 , fh a = 3 x ICC 6 respectively. 


Summary. In this work we have considered the impact of spatially varying stochastic mag¬ 
netic fields and resonant photon-axion mixing on CMB /x distortion. The contribution of mag¬ 
netic field to n distortion depends on the cut-off scale A b and on the damping scale Xd- On 

4 To be more precise, ADMX collaboration did not find any axion in the mass range 3.3p eV-3.5 fj- eV. 
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Figure 2: Sensitivity plot in the parameter space Bq — \b for PIXIE/PRISM [Q] expected value 
on n ~ 5 x 10~ 8 and for the DFSZ axion model. Values of magnetic field spectral index n and 
axion mass fh a are the same as in Fig. [l] 




(a) (b) 

Figure 3: In (a) the exclusion plot in the parameter space Bq—Xb for COBE limit on /r < 9x 10” 5 
for the KSVZ axion model and n = 2 is shown and in (b) the sensitivity plot in the parameter 
space B — \b for PIXIE/PRISM expected value on n ~ 5 x 10 -8 for the KSVZ axion model 
and n = 2 is shown. Values of the axion mass are same as in Fig. [l] 
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Magnetic field (n=-2.9) 

1 q- 9 --Magnetic field and axion (77T a =4.88* IO -5 eV) 

------ Magnetic field and axion (7fT a =1 *10 -5 eV) 

. Magnetic field and axion (m a =3.5* IO -6 eV) 

•10-10 L_m -—.——— —— - --- -————J 

10 100 1000 10 4 io 5 io 6 

Ae(pc) 

(b) 


Figure 4: Exclusion plot in the parameter space Bq — Xb for COBE limit on n < 9 x 10 -5 . 
In (a) the region above the solid line is excluded with no photon-axion mixing included while 
the regions above dashed and dot dashed are respectively excluded by KSVZ and DFSZ axion 
models only. In (b) exclusion plot for magnetic field spectral index n = —2.9 and DFSZ axion 
model is shown. Values of the axion mass are same as in Fig. [lj 
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Magnetic field (n=-1) 

10 _n ----- Magnetic field and axion (7rT a =4.88*10 -5 eV) 

------ Magnetic field and axion (7rT a =1 * 10 -5 eV) 
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Figure 5: Sensitivity plot in the parameter space Bq — A b for PIXIE/PRISM expected limit on 
fj, ~ 5 x 1CP 8 . In (a) the sensitivity plot for KSVZ axion model and n = —2 is shown while in 
(b) the sensitivity plot for KSVZ axion model and n = — 1 is shown. 
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the other hand axion contribution is scale independent as can be seen from Eq. (14). Taking 
into account axion contribution to y, distortion, in general, one finds weaker limits on the scale 
averaged magnetic field Bq in comparison with no photon-axion mixing included. Our main 
results have been shown as exclusion and sensitivity plots in the Bo — Xb plane where the value 
of the chemical potential has been chosen either equal to the upper limit found by COBE or 
equal to the expected limit of PIXIE/PRISM. 

In this work we have considered only resonant photon-axion mixing on generating a non zero 
chemical potential. In the resonant case, the axion mass is not arbitrary but is connected with 
the y epoch redshift, z^. This constraints the resonant axion mass in the range 2.66 x 10 “ 6 eV < 
fh a < 4.88 X 1(E 5 eV [3]. Axions with masses outside this interval make non resonant transition 
into photons. 

The inferiority of axion contribution to y distortions in comparison with magnetic field 
damping is scale dependent as given by Eq. © and Eq. ( fT8] ) . Obviously this is parameter 
depended and certain number of approximations are in order. If we consider magnetic field 
generated by casual mechanism, n > 2 (12], axions with a mass rh a > 3.5 x ICE 6 eV allowed 
by ADMX [ T3] and DFSZ axion model we find that bigger contribution to y distortion with 


respect to magnetic field damping occurs for cut-off scales, A b <172.5 pc and Xb > 4 x 10 3 pc. 
In the Bo — A b plane axion contribution dominates over magnetic field damping for scale 172.5 
pc < Xb < 4 x 10 3 pc. For scales A b — Xb, Bq weakly depend on the spectral index n. For 
example by using Eq. ( |15[ ) we have that for (i<9x ICE 5 , fh a = 2.66 x ICE 6 eV (lower limit 
of resonant axion mass) and DFSZ axion, the upper limit on scale averaged magnetic field is 
Bo < 8.77 x ICE 8 G for n = 2, Bo < 8.49 x ICE 8 G for n = 3 and Bo < 9.81 X ICE 8 G for n = —1. 
If we had neglected the contribution of resonant photon-axion contribution to y distortion, we 
would get Bq < 1.36 x ICE 8 G for n = 2, Bo < 8 x ICE 9 G for n = 3 and Bo < 3.45 x ICE 8 
G for n = — 1. Therefore we can conclude that for values of the parameters as assumed above, 
resonant photon-axion production gives weaker limits on Bo up to a factor in 10 in comparison 
with no photon-axion mixing included. 

Limits on Bo in the case of KSVZ axion model are in general stronger than those from 
DFSZ axion model (for a given axion mass). For example, if we consider axions with mass 
m a = 4.88 x ICE 5 eV (allowed upper limit) and y fixed, contribution of resonant photon-axion 
production to Eq. ( [I5| ) and Eq. ©, is almost marginal. In the Bo — Xb plane the curve 
corresponding to the KSVZ axion model and fh a = 4.88 x ICE ' 5 eV is almost indistinguishable 
from the curve corresponding to the case with no photon-axion mixing included. In general the 
two curves differs from each other up to a factor 1.5 for 172.5 pc < A# < 4 x 10 3 pc. 

A forecast for the future space missions PIXIE/PRISM has been presented. In Fig. [2] we 
can see the level of sensitivity in the case where photon-axion mixing is not included and in the 
case when it is. In the former case we find that, for scales Xb — 415.5 pc, PIXIE/PRISM will 
probe magnetic fields Bo > ICE 10 G while in the latter case is parameter depended. For axions 
with masses fh a > 2.66 x ICE 6 eV we find Bo > 2 nG for DFSZ axion and Bo > 8.8 x ICE 10 
G for KSVZ axion. In the case of axions with masses m a = 4.88 x ICE 5 eV contribution of 
photon-axion mixing is marginal as can be seen from Fig. [2] 

It is important to stress two things. First in this work we did not consider any limit from 
CMB y distortion since at the moment we are currently working on an effective approach of y 
distortion in the case of photon-axion mixing. Second, contribution of magnetic field damping 
to y distortion has been derived by using linearised magnetohydrodynamic equations in the 
photon diffusion limit[5]. Moreover, we have extrapolated both Eq. (15) and Eq. (16) for 
Xb — > Xb ■ A detailed behaviour of Bq around A o requires numerical integration that is beyond 
the scope of this paper but we expect that extrapolated results are very accurate for A b —> Xb- 
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